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Abstract 

In homogeneous universes the propagation of quantum fields gives rise to pair creation 
of quanta with opposite momenta. When computing expectation values of local operators, 
the correlations between these quanta are averaged out and no space-time structure is ob- 
tained. In this article, by an appropriate use of wave packets, we reveal the space-time 
structure of these correlations. We show that every pair emerges from vacuum configura- 
tions which are torn apart so as to give rise to two semi-classical currents: that carried by 
the particle and that of its 'partner'. During inflation the partner's current lives behind the 
Hubble horizon centered around the particle. Hence any measurement performed within a 
Hubble patch would correspond to an uncorrelated density matrix, as for Hawking radia- 
tion. However, when inflation stops, the Hubble radius grows and eventually encompasses 
the partner. When this is realized the coherence is recovered within a patch. In this 
paper, we focus on the case of a massive field with rare pair creation events. However, 
our analysis also applies to cases leading to arbitrary high occupation numbers. Hence it 
could be applied to primordial gravitational waves and primordial density fluctuations. 
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I. INTRODUCTION 

In inflationary models, the large scale structure of the universe is of quantum origin: it 
results from pair creation processes induced by the expansion rate [1, 2]. These primordial 
fluctuations are at the same time very coherent and incoherent, depending on which prop- 
erties one is looking at. On one hand, since the background space is homogeneous and 
since the initial state is vacuum, quanta characterized by different conformal momenta k 
are incoherent. On the other hand, since these quanta are created by pairs, there exist 
(EPR) correlations amongst the two partners in each pair. In fact, at the end of inflation, 
the two-mode states made of k and — k are highly squeezed, so squeezed that one gener- 
ally abandons the quantum settings and works with a classical description of primordial 
fluctuations. Nevertheless, this classical description (in terms of random fluctuations) still 
incorporates the correlations between the k and the — k modes [3, 4]. 

In this article, we analyze an aspect which so far has received little attention, that of 
the space-time distribution of these correlations. This problem is not specific to primordial 
density fluctuations, it also applies to primordial gravitational waves and more generally to 
all pair creation processes in cosmology. To extract the space-time structure encoded in the 
two-modes states it is necessary to introduce wave packets. Indeed, since the background 
is homogeneous, expectation values of local operators are translation invariant and exhibit 
no specific space-time structure. 

The first question to address thus concerns the use of wave-packets for studying pair 
creation. This question has been already addressed to study two related phenomena. It 
was introduced in [5, 6] to reveal the correlations between charged quanta produced in a 
constant electric field and was then applied to Hawking quanta in the context of black 
hole evaporation [6, 7]. In both cases, the space-time structure of the current carried by 
a single pair was obtained. The new ingredient which allows to consider this current is 
the projector II which filters out the final configurations so as to isolate those associated 
with the pair under consideration. Physically, this filtering procedure can be conceived 
as resulting from a local measurement performed at late time (in the case of electro- 
production, it could correspond to detect an electron, or some electric current). One then 
compute {HJ(t, x)) of some local operator J, typically a current. This expectation value 
gives the value of J which is conditional to find the final state associated with II. In 
other words, (UJ(t,x)) determines the space-time distribution of the configurations which 
are correlated to this final state. This distribution shows that to every produced particle 
(described by a local wave-packet) corresponds a well-localized partner which lives on the 
other side of the horizon. In the black hole case, this horizon coincides with the event 
horizon, whereas in the electric case, it is defined by the uniformly accelerated trajectory 
followed by the particle once it is on-shell. 

When applying this procedure to cosmology, we find a similar picture. We shall proceed 
in two steps. We first consider the simple case in which pair creation events are rare. This 
is realized by considering a field whose mass m is much larger than the Hubble parameter 
H. In this case the mean occupation number is exponentially suppressed, of the order 
of e ~ 2 ™/ H . Secondly, we consider the opposite regime wherein the occupation number 
is very large. This is similar to primordial gravitational waves and fluctuations of the 
inflaton field. We shall see that the parameters which fix the space-time distribution of 
the correlations are not related to the mean occupation number. (The latter is determined 
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by the norm of the coefficients of the Bogoliubov transformation whereas the space-time 
distribution of the correlations are fixed by the phase of these coefficients.) Hence our 
analysis applies to arbitrarily high occupation number. In fact it equally applies to two- 
modes states which are described by a classical probability distribution. 

We find that every produced pair of particles results from a dipole of vacuum configu- 
rations which is torn apart by the expansion of the universe and which gives rise to two 
on-shell currents. The creation occurs when the decreasing (physical) momentum is of the 
same order as the mass of the field. The temporal and spatial extension of this region is 
of the order of the Hubble radius H -1 , independently of the mass of the particle. (This is 
reminiscent to pair production in a constant electric field wherein the creation region [6] is 
of the order of the acceleration -1 .) We find that the current carried by the partner of every 
particle lies always on the other side of the Hubble radius centered around that particle. 
Hence any measurement performed within a patch would appear to result from a den- 
sity matrix, as for Hawking radiation when one performs measurements outside the event 
horizon. However, when considering the evolution from a de Sitter period to an adiabatic 
one, the partners will progressively enter the Hubble radius. When they have entered the 
enlarged Hubble radius, local measurements will be sensitive to quantum correlations and 
interference patterns can be obtained. 

In the regime of large mass, the role of the projector n which isolates a given pair and 
the physical interpretation of the corresponding expectation value (HJ(t,x)) are identical 
to those of electro-production: when a heavy particle is detected in a given space-time 
region, one can evaluate the value of the current (or the stress-energy) carried by this 
particle and its partner. In the other regime, when facing the production of a macroscopic 
number of light (or massless) quanta, the procedure should be slightly modified. The 
appropriate way to filter out the final configurations is now achieved by a projector n 
which characterizes a local energy density fluctuation and no longer the particle content 
of the state. Having defined this projector, (UJ(t,x)) tells us the location of the 'partner' 
energy fluctuation. We shall see that the macroscopic character of the occupation number 
does not erase the local character of these correlations. We can thus envisage to apply 
our analysis to physical cosmology wherein primordial gravitational waves and density 
fluctuations are described by highly squeezed two-modes states. This program is only 
sketched at the end of this paper and shall be further developed in a forthcoming work. 

In section 2, we recall the basic steps involved in analyzing pair creation in cosmology. 
We pay special attention to identify the phases which are responsible for the space-time 
structure we are seeking. Then we proceed by an exact treatment of the modes and we 
conclude by an adiabatic treatment which allows to interpret the various results in semi- 
classical terms. In section 3, we show how to introduce wave-packets. We first present the 
results based on the exact solutions by several figures. We then return to the adiabatic 
treatment to explain the origin of their properties. In section 4, we briefly comment on 
the relations between our analysis and the correlations which lead to acoustic peaks in the 
CMB. 
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II. PAIR CREATION IN COSMOLOGY 
A. General formalism and physical parameters 

In this section we provide our notations and discuss the physical meaning of the three 
parameters which govern Bogoliubov coefficients in cosmology. 

Throughout the paper, we consider field propagation in flat Roberson- Walker space- 
times. Their line element is given by 

ds 2 = -dt 2 + a{t) 2 dx 2 , (1) 

with x the comoving position, t the cosmological time, and a the scale factor. The Hubble 
parameter is, as usual, given by H = dta/a. Let $ be a massive complex scalar field 
minimally coupled to gravity. In terms of the proper time, the Klein-Gordon equation 
reads 

d?$ + 3Hd t <i>- -^V 2 $ + m 2 $ = 0. (2) 

To eliminate the first order derivative term multiplied by H, we work with the rescaled 
field 4> = a 3 / 2 ^. In addition, we decompose it into modes of given conformal wave-vector 
k: 

<Kt,x) = J d 5 kj^Mt). (3) 

Then Eq. (2) reduces to a set of decoupled equations 

(d 2 + n 2 )Mt) = o, (4) 

where k = \/kk and where the time-dependent frequency is given by 



«l(t) = . 2 {t) - (|(^) 2 + \\ with c 2 (t) = m 2 + ^ . (5) 



In second quantization, the Fourier components </>k of the rescaled field operator are de- 
composed as 

Mt) = Mt)^ + <PUt)bl k , (6) 

where and fc^ k are respectively the annihilation operator of a particle and the creation 
operator of the corresponding antiparticle with opposite momentum. Moreover, since the 
situation is isotropic, we can and shall work with time dependent modes which 
depend only on the norm of k. To define the creation and destruction operators, one 
needs to choose positive and negative norm solutions of Eq. (2). In terms of the rescaled 
field, the Klein-Gordon scalar product reads 



= i 3 (k-k') (7) 
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where ( , ) is the Wronskian. In conformity with the particle interpretation, we work with 
modes cf) k of unit Wronskians: ((f) k ,(f) k ) = 1 and = — 1- The (time independent) 

creation and destruction operators are defined by 

Ok = ($k, 4>kJ , = - 4>kJ ■ (8) 

In virtue of the equal-time commutation relations 

[k(t)Mlv(t)]=i8 3 Q*-*), (9) 

they obey 

[a k ,4,] = [6k, 6^,] = 5 3 (k-k'). (10) 

Using these operators, the k-th component of the vacuum is defined by the 'two-modes' 
ground state 

|0k) = |0k,p)®|0_ k ,o), (11a) 
a k |0 k ,p) = 0, 6- k |0_ k ,a) =0. (lib) 

In a static universe, this would be the end of the story. However since the frequency 
0, k is time dependent, the notion of positive and negative frequency modes is ambiguous. 
Hence, to set in vacuum at early time Tj n and to read out the particle content at late 
time T out , we must introduce two sets of positive frequency modes of Eq. (2): in modes 
defined at T in and out modes at T out . Because of the homogeneity of the background 
space, the Bogoliubov transformation which relates them reduces to a set of coefficients 
which depend on k only: 

4>T t (t) = a k 4 n (t)+(3* k 4 n *(t). (12) 
By definition, a k and (3 k are given by the (conserved) Wronskians 

o* = (C,C). /$ = -(C*,C) • (13) 

Then, as in Eq. (8), in and out operators are defined by in and out modes respectively. 
These two sets of operators are also related by a Bogoliubov transformation, and they 
define the in and out vacua, |0, in), |0, out) as the tensorial product over k of two-modes 
ground states |0 k ,m) and |0 k ,out) defined as in Eq. (11). 

The normalization of the current (4> k , </>&) = 1 for both in and out modes implies that 

\a k \ 2 - \(3 k \ 2 = 1 . (14) 

This relation reduces to three the number of independent parameters. A convenient way 
to parameterize a k and (5 k is provided by 

a k = e Wk chr fc , (15a) 
p k = e *(»k+2ifc) shrfc . (15b) 



The relationship with the conventional way to describe squeezed states is provided in 
Appendix B. To physically interpret the angles 9 k and ip k requires some care as they do 
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not possess invariant meaning. Instead, r k is easily identified since the mean occupation 
number of out quanta in the in vacuum is given by 

{0in\a° k ut ^a o k f\0in} =5 3 {k- k') \(3 k \ 2 = 5 3 (k- k') sh 2 r fc . (16) 

A similar interpretation of Ok and V'fc cannot be reached. Indeed, by an appropriate choice 
of the arbitrary phases of in and out modes, both k and ip k can always be gauged away 
to zero. In spite of this, as we shall see in the next sections, Ok and 4>k contain physical 
information, given the two sets of in and out modes. To understand this, it is appropriate 
to work with modes such that the phase of 4> k n is zero at Tj n and that (j)° k ut at T out . In 
these settings, the roles of Ok and ipk are easily identified. 

The phase of a k governs the evolution of the modes from Tj n to T out . Indeed, in the 
adiabatic limit, see Section 2.3 for details, one obtains 

k = I dt uj k (t) , (17) 

which is the classical action from Tj n to T out . Thus when considering as initial state 
a superposition of in states with different occupation numbers, and when re-expressing 
this superposition in terms of out states, the changes of the relative phases between the 
various components of the state are governed by Ok- From this fact we learn that when 
the (Heisenberg) state is the in vacuum, k is no longer accessible. This is why it drops 
out in inflationary models when the state of the primordial fluctuation modes is taken 
to be the (Bunch-Davies) vacuum [2, 8]. However when the initial state contains some 
superposition of excited states [9-11], or two inflationary phases [12], the power spectrum 
possesses rapid oscillations which are governed by Ok- 

The second phase ip k contains the information which specifies when the k, — k pair is 
created. To show this is one of the main purpose of Section 2.3. One can already notice 
that ipk is fixed in inflationary scenarios. However this is some how hidden by the (simple 
and appropriate) convention which consists in putting to zero the decaying mode, i.e. of 
keeping only the constant mode before re-entry, when the wavelength of the mode is still 
larger than the Hubble radius [13]. 

In the following, we will be needing the expression of the in vacuum as a superposition 
of out states: 

| 0k , in) = ^- exp f-^a^b^) |0 k , out) . (18) 

|«fc| V a k J 

The derivation of this equation can be found in Appendix B. From Eq. (18) we see that only 
tp k appears in the relative phases between out states with different occupation numbers. 



B. The exact modes 

In a flat de Sitter space, the conformal factor obeys a(t) = jje Ht . In this case, the 
exact solutions of Eq. (4) are given by Bessel functions. The in-modes of unit positive 
Wronskian are given by 

ffW^e-/ 2 ^^), (19) 
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where v = \Jm? / 1 H 2 — 9/4 and where is a Hankel function of the first kind. The 
identification of this solution with an in mode follows from the early time development of 
Hj^(z). For \z\ — ► oo at fixed v [14] one has 
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HgW^J—e^-V e -/2. (20) 



In terms of the conformal time r\ = —e~ Ht , one verifies that <ft£ describes excitations of 
the Bunch-Davis vacuum [8] as it obeys 1 j\fa — ► e~ l ^ v /y/2k. 
The out-modes of unit positive Wronskian are given by 



4>T\t)= -02= M ke~ Bt ). (21) 
\J 2smh(7rz/) 

The identification follows from the asymptotic late time behavior of Ji u . When \z\ — > 0, 
one gets 

/ 1 \i" 

<22> 

Hence, in the late future and up to a constant phase, cj)™* tends to e~ tmt /\/2m which 
corresponds to a comoving massive particle. More details about the asymptotic limits and 
the notion of particule can be found in subsection C. Using |T(1 + iu)\ 2 = 7r^/sinh7r^, one 
verifies that has a unit Wronskian for v > 0. In this paper we shall restrict ourselves 
to v > 0. Hence we shall not analyze minimally coupled massless fields which correspond 
to iv = 3/2. We shall return to this interesting case in a forthcoming work. 

Having identify the modes, the Bogoliubov coefficients defined in Eq. (13) are immedi- 
ately given by the relation 

e _ /2 a, _ s-ftJM zf^M . (23) 

lv sinh(z/7r) 
One gets real fc-independent coefficients given by 

a = ^=J_, = ae~™ . (24) 



VT- e- 



The /c-independence expresses the stationarity of the process and stems from the fact that 
H is constant. Hence a change in k can be absorbed by a harmless shift in time. It should 
also be pointed out that we are working in a 'gauge' wherein both 9^ and ipk are zero. 

Before examining the adiabatic limit, we briefly discuss the space-time behaviour of 
these modes. As such both in and out modes are completely delocalized. However one 
easily obtains the space time behaviour they encode by forming wave-packets: 

/oo 
dk f(k;k) Mt) ^ k{x ~~ x) ■ (25) 
-oo 

We work for simplicity in 1 + 1 dimension and we designate by k the mean conformal 
momentum. Given the phase conventions we have adopted, x corresponds to the position 
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of the 'center of mass' of the pair of quanta involved in the waves packets. This will 
become clear in the sequel. In order to get simple and analytic expressions for </>, we 
choose f{k;k) = 0(k)e~ k / k . This enables us to use Eq. 6.621-4 of [15] 

[ ike -^ m ^{^^j\ (26) 

_ <— > _ 

In Figures 1 and 2, we present the current J = i(f>* dt (f> carried by a wave-packet built 
with out modes in the adiabatic regime, \(3\ <C 1, and for the large occupation number 
\(3\ 3> 1 respectively. The corresponding wave-packets built with in modes exhibit less 
clear features because, for large t, the conformal distance between the particle and its 
partner becomes comparable to the spread of the wave-packet, see Appendix A. For this 
reason, we have plotted the norm of cf> (Fig. 3) in order to show the interferences between 
the two outgoing components. 

In both regimes, at early times, one clearly sees the classical trajectory of the particle 
and that of its partner. These two trajectories are symmetrically distributed with respect 
to x, put to zero in all figures. It should be stressed that wave packets of out modes (or 
in modes) characterize both the semi-classical trajectory of the particle and that of its 
partner. Indeed this pairing in no way depends on the choice of the function f(k, k) but 
results from the creation process itself. It should also be stressed that, as such, the current 
J or the norm |c/>| 2 have no physical meaning. They do not emerge from expectation values 
of operators and should thus be conceived as providing only a pictorial understanding of 
the local character of wave packets. On the contrary, the (more complicate) expressions 
we shall use in Section 3 do follow from expectation values of operators and possess a well 
defined physical meaning. 




FIG. 1: The current carried by a wave-packet built with out modes in the adiabatic regime, for 
\(3\ 2 = e~ 2lx , v — 1, as a function of the conformal coordinate x and the cosmological time t. One 
clearly sees the incoming particle coming from negative x and its partner coming symmetrically 
from positive x. Both follow classical trajectories characterized by their mean momenta k = ±1 
respectively; see App. A. The ratio of the incoming currents is given by — |a//3| 2 . For large positive 
t one just has the unit current carried by the outgoing particle. 
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FIG. 2: The out current as a function of x and t, for k = 1,5 = and v = 0.01, i.e. very far from 
the adiabatic regime. Nevertheless the semi-classical properties of wave packets are robust. The 
incoming currents of respective values |a| 2 and — 1/3| 2 are now of the same order. For large positive 
t one still has the unit current of the particle. 
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FIG. 3: The norm of a wave-packet mode built with in modes for v = 1, k = 1. In addition to the 
incoming particle which enters from negative values of x, there are oscillations for positive x and 
from t ~ 0, when the particle stops to be relativistic. These oscillations arise from the interferences 
between the two out wave packets in the norm | a<t> out + (3<j)* ut \ . (They have been amplified by using 
a non-linear vertical scale.) Deeper in the adiabatic regime, (3 ~ e~^ v <C 1, they are exponentially 
suppressed, whereas in the opposite limit fi/a — * 1 the picture becomes more symmetrical in x, —x 
since both terms in the norm become of the same order. 

C. The adiabatic limit 

The usefulness of considering the adiabatic regime is that it provides a neat semi- 
classical description of the creation process. We shall see indeed that the value of the 
Bogoliubov coefficient [3 is given, as in a tunneling effect, in terms of the exponential of 
the Hamilton-Jacobi action evaluated along a well-defined trajectory. Moreover we shall 
learn from this expression when the creation of the pair occurs, thereby explaining the 
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properties of the former figures. From this identification, we shall also be able to identify 
the modes which are subject to pair creation (and those which are not) when the inflation 
lasts a finite time. 

The adiabatic regime is obtained when the mass of the field m is much larger that 
the Hubble constant H. In this limit, the frequency Qk of Eq. (4) can be replaced by 
the 'classical' frequency uj k = \Jk 2 /a? + m 2 since the differences scale as (H/m) 2 . In 
this regime, the WKB approximation is good. Indeed, its validity is controlled by the 
dimensionless function 



W) = -^ = --A ^~ r - 7T777T I • (27) 



dttok H ( ma 1 
uj 2 m \y k (i +m 2 a 2/ fc 2)| 

The maximum of |A| is given by (i?/m)\/l2/9. It is reached when k/a = my/2, i.e., when 
the physical momentum is of the same order as the mass of the field. These two results 
are exact when H is a constant but are also valid provided H/H 2 <C 1. In the rest of this 
Section, we shall consider only the de Sitter case. The generalisation of the results, such 
as Eq. (45), is obtained by replacing H by H(k) = H(t tp ), where U p (k) is the solution of 
k/a(t) = m. 

It is interesting to notice that X k — > in both asymptotic regimes (t — > ±oo) for different 
reasons: For t — > — oo, the particle goes on the light cone and the mass contribution to 
its energy becomes insignificant. Hence it becomes asymptotically decoupled from the 
expansion because of conformal invariance. In the other regime instead, for t — ► oo, one 
obtains a heavy particle at rest with respect to the cosmological frame and asymptotically 
insensitive to the expansion rate. 

In brief, as in the former subsection, we see that one reaches adiabaticity for asymptotic 
times, for all values of v > to be precise. Hence the identification of asymptotic in and 
out modes is unambiguous. In the adiabatic regime, when v S> 1, one has A&(i) <C 1 for 
all times since A < v~ l . From now on we work in this regime. Then the WKB mode (f>k{t) 

fat) = -j^e^&<n" u ''W , (28) 

is a good approximation of positive frequency solutions of Eq. (4) . This can be seen from 

idtMt) =(l- i^-) oj k (t)4> k (t) . (29) 



Hence it is appropriate to express exact solutions of Eq. (4) as linear superpositions of 
WKB modes: 

Mt) = c k (t)4>k(t) + d k (t)4>* k (t) . (30) 
Owing to the previous discussion, (f>k becomes 1 an exact solution for t — > — oo. Hence, the 



1 One can evaluate the residual effect engendered by the fact that the limit Tin — > — oo has not been 
taken. One finds [11] that the adiabatic vacuum enforced at T; n differs from the Bunch-Davies vacuum 
by Bogoliubov coefficients pT in ,-oc — A fe (T in ) ~ J ,~ 1 e ~ H ' Ti ™~' t p' where t tp (k) is defined by Eq. (43). 
Similarly, the fact of imposing vacuum at some finite late time defines out modes which differ from the 
asymptotic ones by Bogoliubov coefficients /3T ollt ,oo — A fc (T out ) ~ jy _1 e _2H ( T °"* _t 'p). Since the residual 
effects fall off exponentially fast both for early and late times, the notion of asymptotic quanta is well 
defined in the present case. Hence, pair creation amplitudes are also well defined. 
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exact positive frequency solution (j) k m whose phase vanishes at Tj„ (large and negative) is 
given by Eq. (30) with 

c k (T in ) = l, d k (T m ) = 0. (31) 

Similarly, since 4> k is also an exact solution of Eq. (4) in the asymptotic future, it asymp- 
totically coincides with the exact out mode (f) k owt (t), up to an arbitrary phase. We fix this 
phase by requiring that the phase of 4> k out vanishes for t = T out . Hence, for large positive 
t, one has 



Tou t{t) = _J =e if?°<* = j^dt'^it')^^ (32) 



'k 



Moreover, since c(t) and d(t) become constant one can write 

^»(t) = c k (T out ) fa(t) + d k (T out ) fa(t) . (33) 

Together with Eq. (32), this equation gives us the Bogoliubov coefficients between (fij) n 
and (jy[ out . Using their definitions Eq. (13), we get 

a k = 4(T out )e iI ^ tdt '^ t '\ (34a) 
-Pk = dUT out )e- i £r dt '^\ (34b) 

To compute c k (T out ) and d k (T out ), one should determine how they evolve from Tj n to 
T^. Their evolution is completely fixed by requiring that, for all t, one has (see [16] for 
details) 

id t f™ = (<*(<)& - d k (t)4>* k ) ■ (35) 
Then the conservation of the Wronskian gives 

{f k \f k n ) = \c k it)\ 2 -\dkit)\ 2 = ^, (36) 

and Eq. (4) leads to 

Q 9 t UJ i2 /' dt'iU(t') , [ , /0 „ S 

o t c k = —e JT - d k (t) , (37a) 

a rl d tUJ - i2 St dt'w(t') , . , Q „,, 

d t d k = —e JT in 'c k (t). (37b) 

Up to now no approximation has been used in the search of the evolution of c and d. 
However the equations have been recast in a form which is most appropriate to evaluate 
perturbatively non-adiabatic effects. Indeed, in the zeroth order approximation, one ne- 
glects Afe, hence c k = 1 and d k = for all t. This is the usual WKB approximation. To 
first order in the non-adiabaticity, one puts c k = 1 in Eq. (37) and one integrates it to get 

d k {T out ) = [ T ° Ut dt> ^ e -*/£**<*> . (38) 
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Given that A^ <C 1, one can evaluate this integral by a (kind of) saddle point approxi- 
mation. One search for the time t c (k) in the complex plane where cjfc(i c ) = and one 
evaluates the phase of the integrand at that time to get 

d k (T out ) = -ie- 2i %n dtUk(t) . (39) 

The evaluation of the prefactor is rather delicate, see [17] for the details. This adiabatic 
result is valid if the domain [Tj n , T out ] is wide enough so that the saddle time t c (k) is well 
enclosed. If not, one gets boundary contributions which decrease like A/ c (Tj n ) and X k (Ti n ), 
see the former footnotes. Using Eqs. (15b, 34b) we can rewrite dk(T out ) as 

dk(T ou t) = -fte-*" ~ -r k e' 2i ^ +9k) = - e - 2i (^+**) , (40) 

where we have anticipated the fact r k <C 1 in the adiabatic regime. We have also introduced 
the complex phase ip' k = ip k + (V^) hir^. 

In brief, to first order in the non-adiabaticity, and when working with the convention 
that the phase of (p Tin and (p Tout vanish at T, n and T out respectively, 6 k and %// k are given 
by Eq. (17) and 

Yi = -/ dtu k (t)--. (41) 

Jt c (k) 4 

This equation shows that there is a simple and universal expression for M in terms of 
the Hamilton-Jacobi action evaluated along the semi-classical trajectory which originates 
from t c where uj k vanishes. The mean occupation number is fixed by the imaginary part 
of this action. More importantly for us, the real part of ip. (= ip k ) fixes, as we shall see, 
the space-time properties of the correlations induced by pair creation. 

In a flat de Sitter space, all the previous expressions can be evaluated. The location of 
the saddle time given by a(t c ) = H~ 1 e Htc = —ik/m furnishes 2 

<42) 

For later convenience we introduce the 'turning-point' 

1 ffk 

ttp (k) = ^(t c ) = -ln(—) , (43) 
rl m 

where k/a(tt p ) = m. From Eq. (43) we learn that the pair is preferably created when 
the adiabatic parameter A& is near its maximum value, when the particle stops to be 
relativistic. (A pictorial interpretation of the turning point can be reached by adopting 
the somewhat artificial coordinate (t,x) — > (t,y/ax). Then the classical trajectory (see 
Eq. (A7)) possesses a turning point at tt p .) 
Using Eq. ( A4) , we get 

= - [F k (T out ] - F k (t c )} -- = -F k (T out ) (44) 



2 The similarity of this equation with the corresponding saddle point condition governing the Unruh effect 
and black hole evaporation is remarkable, see Eqs. (2.49) and (3.43) in [6]. In all cases, the choice of 
the imaginary sign of the saddle time is such that it leads to exponentially suppressed amplitudes. 
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The real part of ip' k is governed by F k (T out ), the primitive of u>k(t) evaluated at T out . 
This results from the fact that the real part of F k (t) vanishes when evaluated at the 
creation time t c (k). On the other hand, the imaginary part of ip', comes entirely from 
Fk(t c ) = —iTTm/2H which determines the occupation number. In fact one finds 

\d k (T out )\ 2 = \(3 k \ 2 = e^ m / H , (45) 

in agreement with Eq. (24) in the limit m/H 1. 

As a last application of the semi-classical treatment, we show that the identification 
of the creation time t c (k) allows for an evaluation of the decay probability of the vacuum 
when considering the quantization in a space-time box of proper length Lp and of time 
interval AT = Tf — T. The probability that no pair is created is 

\(0out\0in}\ 2 = J]^ = exp(-^ln(l + |/3 fe | 2 )J , (46) 

k V k J 

where one should sum only over relevant modes, that is, modes which contribute to the 
decay of the vacuum given the space-time box. When the size of the box is large enough, 
i.e. Lp ^> 1/H and AT 3> 1/H, one can neglect the edge effects and separate the modes 
according to whether or not their saddle time t c (k) falls within the interval AT = Tj — T, 
see [6] for the same analysis applied to the Schwinger effect. 

The modes with their turning point which falls outside the interval [T,Tj] do not 
contribute since the creation process occurs around tt p (k). Instead, modes with their 
turning point that lies within that interval participate to the above sum. Their number is 

V - 2 [ k(Tf) dk Lp - 2 mLp / fc(T/) dk - 2 mHLpAT (47) 

The factor of 2 arises from the fact that we are dealing with a complex field. The integrand 
of the first integral (= Lp/2ira) is the time dependent density of conformal modes in a 
box of fixed proper length. We then reexpress it as a function of k using t tp {k). The upper 
and lower values of the integral are given by k(t), the inverse function of t tp (k), evaluated 
at T f and T L . 

As one might have expected, the decay of the vacuum is governed by an extensive 
quantity, namely the number of 'relevant' modes in a box of space time volume LpAT. 
When inserting Eq. (46) in Eq. (47), one obtains the equivalent of the famous Schwinger 
formula [6]: 

\(0out\0in)\ 2 = exp Ml + e"™ 2 ^ , (48) 

where E = ma is the constant electric force acting on charged quanta. Notice however 
that the way Eq. (47) is obtained, i.e. through an appeal to exponentially growing wave 
vectors (corresponding to initial trans-Planckian physical momenta), is identical to that 
governing vacuum instability giving rise to black hole radiation [6]. 
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III. SPACE-TIME CORRELATIONS 
A. The absence of structure in expectation values 

When working in the in vacuum, expectation values of local operators such as that 
of the current (J M (x)) or the stress-energy tensor (T^ u {x)) do not depend on x because 
of space homogeneity. Two-point functions such as (J ll (x)J u (y)) or (T^ v {x)T il i l/ i{y)') are 
already more interesting as they correlate field configurations at different locations. Hence 
they are sensitive to the entanglement in the two modes states k, — k. However since all k 
participate with the same weight in these vacuum expectation values, very little structure 
is finally obtained. 

It is a worth exercise to compute this residual structure. Indeed it brings out the 
physically relevant phases irrespectively of the adopted conventions. Moreover it provides 
a neat introduction to the forthcoming analysis of wave packets. To determine this residual 
structure it is sufficient to analyze the Green function as the above mentioned two-point 
functions can be derived from it. For further simplicity we work at equal time and focus 
on its spatial dependence. In the in vacuum Green function is given by 



<O m |0(t,x)^,O)t|O m > = |^A e ^|C(t)| 2 

' ' kx [(Kl 2 + \Pk\ 2 M° k ut (t)\ 2 ~ 23? {aiA(C 4 W) 2 »9.) 



(2vr) 3 

Using the phases introduced in Eq. (15) and the fact that in de Sitter space the norm of 
Bogoliubov coefficients are /c-independent, one has 

d 3 k 



(50) 



(O m |^,x)0(t,O)t| Om ) = |a| 2 |^A e ^ [(l + |^| 2 )|02 ut | 2 -2|^|»{(e^2 Bt ) 2 } 

Two interesting limits can be considered. That governing the adiabatic limit and the 
opposite case in which (3 /a — > 1. We shall restrict ourselves to the first case in this 
section. In this case, the first term in Eq. (50) reduces to the flat result, i.e. the Yukawa 
potential. In 1+1 dimensions, for large mx p , it is given by e~ mXp / ^/mx p where x p = \a(t)x\ 
is the proper distance at time t. In the adiabatic regime, using Eq. (40) , the second term 
is given by the real part of 

_y* e - (e ^ r(( )) 2 = ./^^ye- 2 ^) *»>«>. (51) 

We see that the arbitrary phases associated with Tj„ and T ou t drop out from the product 
e l ^k(f)° ut (t^. This is as it should be since the phase of this product is measurable. It 
determines indeed the sign of the corrections to the vacuum term. To evaluate the integral, 
we use again a saddle-point approximation [18]. One thus searches for the solution of 

8 k ( kx - 2 / dt' co k (t') ) \ k=k * = , (52) 



% ( kx-2 ! 

\ Jtc 



=(*) 

in complex A;-plane. One gets 

x = 2x k (t)\ k=k * . (53) 
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That is, twice the classical displacement introduced in Appendix A, see Eq. (A8). The 
saddle-point approximation is valid well outside the Hubble radius, i.e. a(t)x S> 1/H. In 
this case, k* is always real. It gives the value of k such that, when the particle is found 
at x = at t, its partner is at x at the same time. Thus the function 2xk(t) gives the 
(mean) distance between the particle and its partner when they are created from vacuum. 
This will be clarified in the next Section. When xa(t) » 1/H, one asymptotically gets 
k*x ~ 2m/ H. Then the second term of the Green function behaves as 



Contrary to the Yukawa term, the amplitude of this term decreases only as the inverse 
of the conformal distance x. (This slow decrease arises from ^A^, the width around the 
saddle point, given in Eq. (A10).) Therefore, the term linear in (3 dominates the Yukawa 
term in Eq. (49) when the proper distance is larger than ir/H. 

In brief, in the adiabatic limit, the two point function contains the usual Yukawa term 
and a small additional term which is due to pair creation processes. The latter does not 
decrease like an exponential because, for arbitrary large spatial separation, there is always 
a pair of onshell quanta which connects the two points since Xk(t) diverges as k — ► 0. From 
Eq. (54), we see that the large distance behaviour contains oscillations whose wavelength 
grows logarithmically. This absence of a distinct pattern is due to the fact that all values 
of k contribute to the Green function evaluated in the in vacuum. 

The lesson of this analysis is that this Green function is not the right object to unravel 
the space-time correlations. To obtain these correlations, one must introduce a weight 
in k space so as to keep only a finite range of conformal momenta. This can be done 
by considering wave-packets. To understand how to introduce wave-packets in second 
quantization, it is appropriate to pose a physical question whose answer will provide one. 
Consider the following question: knowing that we start from vacuum and that a particle 
has been detected at some late time to around xq, where should one look for its 'partner' 
? The appropriate formalism to answer this question consists first, in introducing the 
projector II which is associated with the detection of the particle, and second, in computing 
the value of operators which are conditional to this detection. This was first done in [5] 
using the formalism developed by Aharonov et al. [19]. 



When a particle is detected at to,xo with some momentum k, the state of the field 
is 'reduced'. This can be viewed from an axiomatic point of view, or better, a la von 
Neumann, from the interactions between the field and some additional quantum system 
which acts as a particle detector. This has been explained with details in [5, 20] and will 
not be repeated here. In one word, the additional system can be thought to be a bubble 
chamber. Hence the detection of a particle is the recording of a track in the chamber from 
which one can deduce both the location and the mean momentum k of the particle [21]. 

For simplicity we consider here only the case of rare pair creation events. Hence we 
can restrict the analysis to one particle states. We thus assume that a particle is detected 
around some xq, at a large future time to- Since the detection occurs for to S> t tp (k), where 




(54) 



B. The projector 
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t tp is defined in Eq. (43), one must use the out basis to characterize the local excitation 
which triggers the external system acting as a particule detector. In (1 + 1) dimensions, 
the state of the detected particule is thus specified by 

\ps) = J dkf k a o k ut ^\0 p out) , (55) 

where the function gives the probability amplitude to find a particle with momentum 
k. A subscript p (a) added to a ket means that it belongs to the particle (antiparticle) 
sector in Fock space. Notice indeed that in Eq. (55) nothing has been said about the 
anti-particle states. It is therefore convenient to introduce the projector associated with 
the detection of \ps). It is given by 

n ps := \ps)(ps\ (g> I a 

= (J dkf k a° k ut ^\O p out) (O p out\ J dk' f^afft) ® I a , (56) 

where I a is the identity operator on the antiparticle sector in Fock space. 

This projector contains all the information about the detected particle and its partner 
when the Heisenberg state is the in vacuum. To show this let us apply H ps on |0m). Using 
Eq. (18) and developing the exponential to first order, we get 

H ps \0in) = — 7=|ps) <8> / dk' f k , — —b^_ k ,\O a out) = — \=\ps) ® \partner) . (57) 
y Z J ay \J Z 

From Eq. (57) we see that the specification of the particle state \ps) univocally fixes the ket 
\partner), the state of its partner 3 . This is due to the EPR correlations between particle 
and antiparticle out states which are present in |0j ra ). The projector Yl ps gives also the 
probability to find the chosen particle: 

P ps = (0m|n ps |0m) = \(0 in \0 out }\ 2 [ dk\f k \ 2 \^\ 2 , (58) 

J Oik 

which is simply the weighted sum of the probabilities to get a particle with momentum k. 



C. The conditional value of the current 

The second element of the analysis of [5, 7] is the notion of the expectation value of 
an operator which is conditional to the detection of the particle described by the ket \ps). 
Hence instead of considering, as we did in the former Section, the expectation value of 
in the in vacuum, see Eq. (49), we now study its conditional value: 

,. (t \ / / </ /xf\ (0m|n ps (/.(x)(/)(x / )t|0in) 



3 Eq. (57) is exact, i.e. valid for all values of f3/a even though we are presently working in the limit 
f3/a < 1. 
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The denominator is the probability given in Eq. (58). The numerator is given by 



/ 



0k' 



dkdk'f£f k ,(-—)^0out\b o ^,a% Lt (l)(x)4(x'y\0in} . 

Oik' 



Using the Bogoliubov transformation a\ n = ctka° k ut + Pkb°™l , the commutation relation 
[b™ k \ b%] = a* k 5{k - k'), we get 



(60) 



{Qout\b™l,a° k ut (i>{x)<P{x'y\Qin) 



OkOik 

0k 



■« iO (z)C;0O<0out|0in> 



oik 



5(k - k')(0out\ct)(x)(l)(x'^\0in) . (61) 



We have added the subscripts p and a in order to easily identify the wave functions which 
refer to the particle or the anti-particle. (Notice that in the presence of an electric field, 
they would be different functions.) Using Eqs. (60,61), the conditional Green function 
Eq. (59) becomes, 



1 



{(t>{x)(t>{xy) cv = G F (x,x') + — 



ps 



dk^cj) k n ;{x') 



[dk'fk,^- 
J a k ,a k ' 



where 



Gf(x, x') = 



(Oout\(f)(x)(f>(x'^\0in) 
(0oitt|0m) 



(62) 



(63) 



is the in-out propagator properly normalized. 

In brief, the conditional value of the Green function splits into a background term 
which is independent of the chosen particle and a term which is specific to what has been 
detected. 

In the next subsection, we shall use the adiabatic treatment to see how the various 
semi-classical elements determine the properties of conditional values. Before doing so, 
we present by several figures the conditional value of the charged density obtained with 
the exact modes of Section 2. (We could have worked with the conditional value of the 
energy density (= (nT o(i, x))/(II)) but since it is more complicated we proceed with the 
current.) In perfect analogy with Eq. (59), the conditional value of the charge density is 
given by 



{J(t,x)) c 



(Om|n ps 0(t, x)idt (f)(t, x)'\0in) 



(0in\H ps \0in) 

It can be straightforwardly obtained from Eq. (62), and it also splits into two terms 



(64) 



(J(t,x)) c 



(0out\J(t,x)\0in) 1 



(65) 



(0out\0in) 



P, 



ps 



dk^;(x) 

Oik 



i d t 



dk'f k 



a k ,ak' 



The in-out term will not be analyzed here as it presents no structure. We focus instead 
on (J(t,x)) ps which is specific to the selected state \ps). (J(t,x)) ps is complex function. 
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At first sight, this seems awkward. However, after a little thought one understands that 
this must be the case. Indeed, the imaginary part of (J) ps governs the linear change of 
the probability P ps to find the chosen particle when one modifies the electric field [6]. 
(Similarly the imaginary part of the conditional value of T„ v governs the change of P ps 
when one modifies the background geometry, see Eq. (33) in [7].) 

The main features of the conditional value of the current are as follows. First, both the 
real and imaginary part of {J(t, x)) ps vanish for t <C tt p = H^ 1 ln(kH/m), much before the 
moment when the /c-particles stop to be relativistic. In other words the ^-configurations 
are still as in vacuum. Second, {J(t,x)) ps gives rise to two well-localized (unit) currents 
which follow the trajectory of the particle and that of its partner. (The fact that these 
currents are unity follows from the presence of the denominator in conditional values. The 
interpretation is clear: when a particle has been detected, the conditional currents are ±1.) 
Third, these two semi-classical currents emerge from wild oscillations which are confined 
in the creation zone, i.e., in a space-time domain of extension 2/H x 2/H centered around 
the turning point U p (k). In this region the imaginary part of (J(t,x)) cv does not vanish 
but oscillates with the same amplitude as the real part. To explain the origin of these 
properties it is very useful to return to the adiabatic treatment. 



2 




FIG. 4: The real part of the conditional current (J(t,x)) ps for v = 3 as a function of the proper 
coordinate X and the cosmological time t. The Hubble radius and the Hubble time are equal to 
one. We have used a nonlinear scale on the z-axis to tame the wild oscillations in the creation 
zone around the turning point t tp ~ — 1. We clearly see that the current vanishes in the past of 
the creation region and that two semi-classical currents emerge from it. Moreover, they propagate 
along classical trajectories, see Figure 8. 



D. The adiabatic treatment 

It is easy to understand why (J(t,x)) ps vanishes in the past of the turning point. It 
follows from the fact that (J) ps is the product of two different wave packets, see Eq. (65): 
one is associated with the incoming particle and the other to its partner. These wave 
packets are centered along classical trajectories which diverge from each other, see Figures 
1 and 7. Hence their product vanishes asymptotically in the past. 
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FIG. 5: The imaginary part of the conditional current presented in the former Figure. We have 
again used a nonlinear scale on the z-axis to tame the fluctuations which are of the same order 
as the real ones. The imaginary part takes significant values only in the creation region. This 
means that once the particles are on-shell, the conditional current (J(t, x)) ps behaves as a classical 
current. 




FIG. 6: The real value of the conditional current far away from the adiabatic regime, for v = 0.01, 
as a function of X and t. Even though the wave packets are much broader (this is hardly surprising 
since the mass is now much smaller that the Hubble radius = 1) one sees that the properties of 
the conditional value are robust. It still vanishes in the past of the turning point and still gives 
rise to two semi-classical currents following classical trajectories. 

To understand the properties in the future of the turning point, namely why it describes 
two localized and real currents, it is appropriate to work with out modes since they are 
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WKB in that period. Using Eq. (12), the numerator of {J(t,x)) ps is 

i d t 



oik 



\(*k\ 



J dkf k <fi£{x) 



i d t 



[ dkf k \^\ 2 ^(x) 

J OL k 



+ 



dkf k ^% a (x) 



i d t 



! 
I 



dkfXf(x) 



i d t 



dkf^out {x) 
Ok 



a. 



dkfk^Ux) 



OLk 



i d t 



I 



dkf k \^\ 2 4>Z(x) 



Oik 



In Eq. (66a), in the first bracket one finds the wave function of the particle, 

F p (t,x;t ,x ,k) = {0 p out\4>{t,x)\ps) = J dk f k e ikx <j) T k ° ut {t) . 
To evaluate it, we use for simplicity a gaussian wave packet 



(fc-fc) 



fk = e I*' 2 e 



y- p -ikx 0p ~i.u° ut dt w *w 



(66a) 
(66b) 
(66c) 
(66d) 

(67) 

(68) 



The phases are chosen so that the wave- function is maximum at (to, xq). We proceed again 
with a saddle-point analysis. Up to an irrelevant phase, we get 



F p (x,t) 



2vr 



d 2 k $(k* 



where 



$(fc) = + ik(x - x ) 



The saddle-point k* is the solution of 

k* -k 



dk lk 



2a 2 



+ i(x - xq) 



t 

i I dt' u k (t') . 
ho 



[ dt' d k u k {t') 

J t 



(69) 



(70) 



(71) 



If one neglects the spread, i.e. a — ► oo, k = k, one recovers the classical equation of motion 
Eq. (A7). The role of the spread is to damped the wave packet when k* differs from the 
mean momentum k. We thus look for a solution of the type k* = k + 5k. Expanding 
Eq. (71) around k to lowest order in 5k, we get 



5k = iTZ\ (x-x - Ax k (t , t)) , 
where Ax|(to,i) is the classical displacement from to to t when k = k. 

-9 k ^k\k = ^+id k Ax k (t ,t)\ k= - k , 



2E|(t) 



2a 



(72) 



(73) 
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gives the time dependent spread of the particle's wave packet. When expanding <&(&*) 
around k one gets 



F P {x,t) ~ \j^^f m ex P (-2S?[x - xo - Ax s (t ,t)] 2 ) • (74) 

As expected, the wave-packet is a gaussian centered along the classical trajectory charac- 
terized by the mean momentum k. 

The second bracket in Eq. (66a) is F p (x,t) multiplied by a k-independent real factor, 
\§^\ 2 - Hence the first line of the r.h.s in Eq. (66) gives the current carried by the wave 
packet described by the asymptotic state \ps). This explains why asymptotically (J) ps is 
real and concentrated along the trajectory associated with the chosen particle state. We 
shall see below that the same reasoning will apply to the partner's wave function. The 
role of the normalization factor \§^\ 2 is to guarantee that the current carried by (J) cv is 
one, thanks for the presence of the denominator P ps in Eq. (62). 

In the third and the fourth lines of Eq. (66), one finds products of the particle and 
the antiparticle wave-packets. For late times they vanish in the adiabatic limit because 
the classical trajectories do not overlap. Their role is to guarantee that (J) ps vanishes in 
the past, and to contribute to the complex oscillations of large amplitude (= \a/(3\) in 
the creation zone. To show that there is no significant overlap in these products, we now 
compute the properties of the wave-packet of the partner which appears in the second line 
in the r.h.s. of Eq. (66). 

In Eq. (66b), we find the wave function of the partner and its complex conjugate. 
Hence the current they describe is real. (In fact, this follows from causality: when t, x 
is separated from the support of n by a space-like distance, n and J commute, hence 
(J) P s is real.) Having understood the asymptotic reality properties of the conditional 
current {J) ps , it is now interesting to determine the trajectory followed the partner and 
what are the parameters which determine it. As for the particle, it is appropriate to use 
a saddle-point approximation. Explicitly, one has: 

F a (t,x) = {O a out\0(t,x^\partner) = f dk f* k ^e~ ikl : (f)f ut (t) = I dk e* (fe) , (75) 

J a k J 

where 

(fr — fr) 2 fto 
^{k;x,t,x ,t ) = 4 2 ~ ik(x - x ) + i2ip' k + i / dt' u k (t') , (76) 

{k-kf I no r t 



4a 2 



-ik(x-x )-i\ f° dt'u k (t') + [ dt' u k (t'U . 

\Jt c (k) Jt c (k) I 



In the first line we have used the complex phase ip' k introduced in Eq. (40) and we have 
identified T ou t with the late detection time to- The unusual sum of actions we obtain from 
ip' k has a clear meaning. It gives the phase accumulated from the detection of the particle 
at time to to the location of the partner at time t when going backwards in time so as to 
pass through the creation process which occurred at t c (k). 

The behaviour of F a is governed by the saddle-point k* defined by ^\k* =0. As in the 
previous analysis, when one searches for the characterization of the wave function close 
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to its maximum, all quantities can be evaluated for k = k. The condition ^-|^ = thus 
gives the trajectory of the partner: 



X a (t) =X + d k 



2^ k + J t0 dt' u k {t') 



\k=k ■ 



(77) 



From this equation, the classical equation of motion of the particle, Eq. (A6) evaluated 
for k* = k, we get 



x a (t) - x p (t) = -2d k / dtu k (t)\ k , 

JtrXk) 



(fc) 

= -2x k (t)\ k=k . (78) 

This is the central equation of this Section. It tells us that, at time t, the (mean) conformal 
separation between the fe-particle and its partner only depends on the HJ action evaluated 
from the creation time t c (k) to t. None of the details of the wave packet enters in it. This 
result agrees we what we found in Eq. (53). Notice that it would have been also obtained 
if we used the conditional value of energy density rather than ( J(t, x)) ps . Moreover, when 
working outside the semi-classical regime, one still reaches the same conclusion, see Figure 
6. 

Expanding the solution around k, one gets 



Fa{x > t] " V WW) em exp (" 2S ^ (x " Xa{t)? ) • (79) 

As in Eq. (74), we obtained the semi-classical wave function evaluated at k = k times 
gaussian centered along the classical trajectory of the partner. They are however some 
differences. First e*( fe ) contains a factor e - nm / H = (3/ a. Hence the current carried by the 
second line in the r.h.s of Eq. (62) is equal to —1 when taking into account the denominator 
P ps . It is thus equal and opposite to that carried by the particle wave function in the first 
line. Second, the interesting novelty concerns the spread. It is given by 

- i2d k x a (t)\ k=k = ^-i [A fc (t) + A k (t )} \ k= - k , (80) 



2T 2 (t) 2a 2 ' "^'^ 2a 2 

k.a x ' 



where A k (t) is given in Eq. (A10). As in Eq. (73), it is the 'susceptibility' of the HJ 
trajectory which determines the time dependence of the spread. The novelty is that this 
spread now depends on the second derivative of tp' k . Hence it is intrinsic to the pair creation 
process. From it we can derive the notion of minimal wave-packets. Consider S k a (to), the 
spread in x of F a (x,t) at time to 3> tt p (k). It is given by 



qjo) ~ -HM ~ i/^ + Ia^o) • (81) 

Its minimum value is obtained for a~ 2 = 2A k (to) ~ 2m/Hk 2 , see Eq. (All). Then 
k 2 5 ka = 8m/ H for all k and for late time. We thus see that the otherwise arbitrary 
spread in k (= a) is now fixed by the asymptotic value of the second derivative of the HJ 
action evaluated from the creation time t c (k). In this respect, minimal wave packet are 
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intrinsic to the creation process. Moreover, since they cover a space-time zone of proper 
extension At x Axp = (mH)' 1 , the number of non overlapping minimal wave packets in 
a space time box correctly gives the number of 'relevant' modes which govern the vacuum 
instability in Eq. (47). Hence they provide a natural basis for describing pair creation 



processes 4 . 



This result also justifies the fact that we neglected the product of wave functions of 
the particle and its partner in Eqs. (66c, 66d). Indeed, the asymptotic separation of the 
classical trajectories of the particle and its partner is x a (to) — x p (to) ~ 2m/ Hk. Thus one 
has 



x a (to) - x p (to) 



(82) 



5k,a V H 

Therefore, for minimal wave packets and in the adiabatic regime m/H S> 1, the asymptotic 
value of the overlaps in Eqs. (66c, 66d) are exponentially small. 

Finally we characterize the size of the space-time region in which the conditional current 
( J) ps does not behaves semi-classically. Its location in time is centered near the turning 
point time Htt p = \n(kH/m) and has a duration of about 4 Hubble times. For larger time 
lapses, the non-adiabaticity of the propagation, governed by Afc(i), decreases exponentially 
fast. Hence non-adiabatic effects are concentrated in that lapse. The spread in space of 
this region is then given by the distance between the two partners at t tp : 



V2_ 

Ha(t tp ) 



Axk,creation=2Xk(t)\ t=ttp ( k ) = 2 — — - . (83) 



Hence the corresponding proper distance is given by Axp = 2y/2/H. Thus, at the creation 
time and for all values of k, the partner is outside the Hubble radius centered around the 
particle by a factor of 2s/2, see Figure 9. As long as the inflation goes on, the two members 
in each pair stay outside the Hubble radius. 

Hence, from all measurements performed within a Hubble patch, one would conclude 
that the created particles are described by an incoherent density matrix, as in the case 
of Hawking radiation when performing measurements far away from the black hole [6, 8]. 
However, the difference with black hole physics is that the Hubble radius grows when 
inflation stops. Hence progressively the particles will join their partner within the same 
Hubble radius. Then, the coherence of the in vacuum is recovered within a patch and 
interference patterns can be obtained. 

IV. LINK WITH PHYSICAL COSMOLOGY 

When applying the preceding analysis of conditional values to primordial gravitational 
waves or primordial density fluctuations, one faces new features. In this Section, we briefly 
present them. A more detailed analysis will be presented in a forthcoming work. 



4 This is again very similar to what is found when studying the Schwinger effect [6]. In fact the whole 
analysis we just performed can be applied as such to study that case. One should simply work in the 
homogeneous gauge At = 0, A x = Et, and replace our frequency uj\ — m 2 + k 2 /a(t) 2 by the relevant 
one: col = m 2 + (k - Et) 2 . 
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The first novelty arises from the fact that primordial waves are described by massless 
minimally coupled fields. Hence, their frequency Qk becomes imaginary when the modes 
exit the Hubble radius. Thus out modes no longer exist. This can be seen from Eq. (22) 
and the fact that v is now imaginary. Instead in modes are still perfectly defined as in 
Eq. (19) because in the past the physical momentum still becomes much smaller than 
the Hubble radius. To obtain well-defined out modes (which are necessary to specify the 
state at late times) one should stop inflation and consider the following adiabatic era so 
as to let the modes re-enter the Hubble radius and start to re-oscillate. This does not 
raise any new conceptual difficulty but the necessity of dealing with two different periods 
of expansion prevents to obtain a simple and analytical description of modes throughout 
their evolution (unless one neglects, as usually done, the decaying mode.) 

The second novelty concerns the occupation number (= |/?fc| 2 ). It becomes extremely 
large during the period wherein the modes do not oscillate. Thus one inevitably works in 
the limit 'au — > 1. In this limit, Eq. (50) becomes 



(O m |0(t,x)0(t,O)t|O m ) = 2L4| 2 -g- \<fi*(t)\ 2 (1 - cos2^(t)) , (84) 

since in inflation one has |afc| 2 = \A\ 2 /k 4 when the slow-roll condition dtH/H 2 <C 1 is 
satisfied [2]. We have introduced the time dependent phase tpk(t) given by 

As in the adiabatic case j3/a <C 1, see Eq. (50), ipk(t) governs the spatial properties of the 
Green function. [Notice that this phase also governs the spherical harmonic coefficients 
Ci(t r , t e ) which characterize the two-point function at emission time t e as seen at reception 
time t r . They are given by [2] 



Ci(t r ,t e ) = 



— — / dkk j l (kArj) — 

9 Jo a(t e y 

\A\ 2 (An) 2 f°°dk , 2 1 - cos2V> fe (t e ) 



a(t e f 



f°°dk 9 „ A . 1 cos2ib k (t e ) , N 

I T jhkAr]) (86) 



where A77 = rj(t r ) — rj(t e ) is the lapse of conformal time and where ji is the spherical Bessel 
function. We have also used the fact that in a radiation dominated universe, the rescaled 
field <fi introduced after Eq. (2) satisfies k\4>^ ut (t e )\ 2 / a(t e ) = 1/2 when m = 0. When the 
function which multiplies jf in the integrand is fc-independent, 1(1 + 1)C; is independent 
of I. Hence, the dependence in I of 1(1 + 1)Q is governed by Y>fc(i e )-] 

The difficulty related to the high occupation number only concerns the use of condi- 
tional values. Indeed, it is now meaningless to specify the occupation number as we did 
it in Eq. (56). One should now work with a new projector which specifies that a classical 
wave has been detected. This projector should thus be built on the field amplitude itself. 
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Instead of Eq. (56), it could be now of the form 5 

IL = Jdk> f k ra\° ut J dk P k d° k ut . (87) 

Perhaps the simplest way to understand its meaning is to consider the interaction of the 
field cj)(t, x) with an additional system which acts as a wave detector, and to work to second 
order in the coupling by letting the dynamics proceed. We refer to Section 3. A of [20] for 
an explicit example. 

Using the new projector, in the place of Eq. (62), one obtains 



(Oin\fl4>(x)4>(xy\Oin) 
(0m|ft|0m) 



in I 



1 

P 



J dk'M-(3* kl )cf>™ ja (x) [| dkr k aM^(x') 



,(88) 



where P is the probability to detect the chosen wave. It is given by P = (0m|ft|0m). 
As in Eq. (62), the conditional value splits into a background term and a term specific 
to the wave which has been detected. We note two differences. The background term is 
now the Green function evaluated in the in vacuum. This results from the fact that the 
new projector specifies much less the final configurations that the first one did. Indeed 
the latter stipulated that only one particle was found, hence the appearance of the out 
vacuum in Eq. (62). The other novelty concerns the second term. Each integrand has 
been multiplied by |«fc| 2 = n& + 1, which results from the Bose statistics of the field. The 
crucial point is that the presence of this new factor contains no phase. Hence the locus of 
constructive interferences is unchanged with respect to we had in Eq. (62). Thus the results 
found in Section 3.4 still apply, no matter how large is the occupation number. In fact the 
analysis would also apply to configurations which are described by a classical probability 
distribution. In this case however, one would obtain conditional values which are real. This 
should cause no surprise since the correspondence between quantum two-point functions 
and classical correlators is obtained when using anti-commutators in quantum settings. In 
our case, the 'classical' conditional two-point function would be given by Eq. (88) when 
replacing Ilcf)^ by the anti-commutator {ft, <jxft}+. This leads to a real conditional value 
which coincides to that one would obtain using classical settings. 

In brief, the usefulness of considering conditional values, and not only expectation 
values, still applies in the classical limit. As in the adiabatic regime, conditional values 
give us the spatial distribution of correlations relative to a restricted set of configurations, 
those which have been selected by the projector. It is this filtering procedure which 
brings much more detailed information as it eliminates the washing out mechanism present 
in expectation values. One can thus envisage to compute the conditional values of the 
harmonic coefficients as they would contain more information as well. They would now 
depend on m, the projection of / on a given axe. They are given by the usual definition 



5 It should be clear to the reader that this is not the only possibility. Since we are working in the classical 
limit, one could equally work with an anti-commutator. When the differences between two options are 
given in terms of commutators, they are irrelevant in the classical limit since they are order 1 whereas 
anti-commutators are order of \f3\ 2 S> 1. 
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when replacing the expectation value of (/>(t e ,x)</>(i e ,x') by its conditional value. The 
relevance of applying these concepts to observational data is presently under consideration. 

V. CONCLUSIONS 

In this article, we have analyzed the space-time distribution of the correlations which 
are induced by pair creation processes in cosmology. Our analysis applies both to the 
adiabatic regime, where the creation rate is exponentially suppressed, and to the opposite 
regime where the mean occupation number can be arbitrary large. The reason is that the 
space-time distribution is determined by products of the form <f>\ out f3k/&k, see Eq. (49) 
and Eq. (84). Irrespectively of the occupation number (= |/3||) it is thus the ^-dependence 
of the phase of (3k / a k which matters. 

However the space-time distribution exhibited in expectation values shows no specific 
properties because all momenta contribute with equal weight in the vacuum. Therefore, to 
isolate the correlations associated with a given pair of particles, we studied the value of the 
current which is conditional to the detection of one member of that pair, see Eq. (65). We 
found that each pair gives rise to two semi-classical currents which are localized along the 
trajectories of the particle and its partner. These currents emerge from wild oscillations 
in a creation zone of space-time extension ~ A/H x A/H which is centered around the 
time at which the particles cease to be relativistic. We then showed that the properties 
of these conditional values are intrinsic to the creation process in that they follow from 
the Hamilton- Jacobi action evaluated from the creation time t c (k). This action and its 
derivatives with respect to k govern indeed the statistical weight of amplitudes, Eq. (39), 
the space-time separation between the particle and its partner, Eq. (78), and the spread 
of their wave functions, Eq. (80). 

Even though this semi-classical analysis is a priori restricted to the adiabatic regime, 
the results it explains are also found very far from it, as shown in Figure 6. In particular 
each pair of modes still gives rise to a pair of localized energy density fluctuations whose 
properties are fixed by the creation process. Thus the space-time information encoded in 
conditional values is still richer than that of expectation values where some details are 
washed out by the averaging procedure. The robustness of the space time properties far 
from the adiabatic regime is an interesting fact which deserves further study. To clarify 
this point, to analyze what happens when dealing with massless minimally coupled fields 
(i.e., the relevant case for describing primordial gravitational waves), and to determine the 
possible use of conditional values in analyzing primordial spectra are three motivations 
for further work. 
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APPENDIX A: THE HAMILTON- JACOBI ACTION 



In this Appendix, we compute the various integrals obtained when working in the 
adiabatic approximation. They are all expressed in terms of the Hamilton-Jacobi action 
and its derivatives. For simplicity we work in a flat de Sitter space in 1 + 1 dimensions. 

In the RW metric Eq. (1), the Hamilton-Jacobi equation reads [22] 



-{d t Sf + ^(d^S) 2 + m 2 = Q 



(Al) 



where k = d x S is the comoving momentum. Since it is conserved, the action separates as 



S(x,t;x ,to) = k - (x - x ) - S k (t,t ) , 



where 



{d t S k f = u k {tf 



+ m 2 . 



The solution with positive energy which vanishes at to can be written as 

rt 



S k (t,t ) = [ dt' u k {t') = F k {t) - F k (t ) , 

J to 



where the primitive is given by 

F k (t) = - 



m 
H 



— \J\ + u 2 — arcsinh(u) 



(A2) 



(A3) 



(A4) 



(A5) 



and where u(t) = ma{t)/k. 

The classical trajectory keeping the end-points fixed is given by 

dkS^tjxo^o) = 0. 
It gives x(t) = xo + Ax c /(t,io)- The classical displacement from to to t is 

sgn(k) 



(A6) 



Ax d (t,t ) = £^3 k u, fc (t')= 



m 2 a(ty 
k 2 



= x k (t) - x fc (to) • (A7) 



J t 



In preparation for the semi-classical description of pair production, we have introduced 

sgn(k) 



x k {t) = d^F k {t) 



Ha{t) 



1 + 



m 2 a(t) 2 
k 2 



(A8) 



As we shall see, it will determine the position of the particle in a pair with respect to the 
'center of mass' of that pair, called x in Section 2.2. In other words 2x k (t) is the conformal 
distance between the two members in all pairs of conformal momenta k, see Section 3.4 
and Figures 7 and 8. Its asymptotic behaviour is 



for t 



-oo, x k ~ 



sgnjk) 
Ha{t) 



-sgn(k)oo 



for t — ► +oo, x k ~ — sgn(k) 



m 1 
Hk 



(A9a) 
(A9b) 
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Finally, wave-packets tend to spread because of the non-linear dependence of Sk with 
respect to k. In the semi-classical approximation, the growth of the spread is governed by 
the susceptibility of the classical trajectory with respect to k: 



m 2 i 



Ak{t) = ^ Xk = m (ky a * + m*)v* • (A10) 



Asymptotically, one finds 



for t — ► — oo , Afc — > , 

in 1 

fort^+oo, & k ^——. (All) 




FIG. 7: The continuous line corresponds to the Hubble radius centered around x = 0. x is the 
conformal distance and t the cosmological time. In dotted and dashed lines, we have drawn two 
pairs of classical trajectories with respective comoving momenta k — ±1 and k — ±l/e. In each 
pair, the two trajectories have been placed symmetrically with respect to x = and the arbitrary 
distance between them has been fixed by the result delivered by the pair creation process, that is 
2x k {t) of Eq. (A8). 



APPENDIX B: SQUEEZING FORMALISM 

In this appendix, we make contact between the Bogoliubov formalism and the squeezing 
formalism. 

The in and out annihilation operators and vacua are respectively related by 

< = Sa^S* , (Bla) 
|0,m) = S\0out) , (Bib) 

where S is the scattering-matrix is defined by 

S = U(T m ,T out ) . (B2) 
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FIG. 8: The same trajectories as in the former figure but now represented in terms of the proper 
distance X = a{t)x measured from x = 0. In our units, the Hubble radius is equal to one. 
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FIG. 9: The trajectory corresponding to k — 1 now represented in terms of the proper distance 
from the asymptotic position of the particle with positive momentum. The Hubble radius is still 
equal to one. In the remote past, the partner's trajectory hugs the outer side of the horizon, as in 
black hole radiation [6]. 



U(t) is the evolution operator, the time ordered exponential of — i fJp_ dt H(t). 

In the squeezing formalism, the evolution operator is decomposed into the product of 
a rotation operator R and a squeezing operator S, see [23]. Since the components of the 
field decouple, one can thus work at fixed k and write 



(B3) 



When expressing these two operators in terms of out creation and destruction operators, 
one has 



J R k (0)=exp|-^(aj { a k + 6 t _ k 6_ k )} , 
S k (r, 0) = exp {r (e"^6_ k a k - cc) } . 



(B4a) 
(B4b) 
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Their action on annihilation and creation (out) operators are: 

R{9)a k R^(9) = e*V, R(9)b_ k R,\9) = e* e 6„ k (B5) 

and 

S{r, (f>)a k S\r, <j>) = ch(r)a k + e i2 *sh(r)6 t _ k (B6) 
We now give the expression of the in vacuum in term of out states. Starting from 

|0 k , in) = R k (9 k )S k (r k , fa) |0 k , out) , (B7) 
and using the general relation 

R\0)S(r, <f>)R{9) = S(r, + 9) , (B8) 

one obtains 

|0 k , in) = S k (r k , (p k - 9 k ) |0 k , out) . (B9) 
One can further decompose the squeezing operator into the product of three operators: 

SUnAk) = ^exp{-e +l2 ^thr fe aj^} 

xexp(-21n(chr fe )(aj c a k + 6^ k 6_ k )) exp( e - i2 ^a k 6_ k ) . (BIO) 
Regrouping Eq. (B9) and Eq. (BIO), one finally gets 6 

|0 k , in) = JL exp (- e +i2 ^- e *hhr k 4&L k ) |0 k , out) (Bll) 

To make contact with Eq. (18), we first identify the Bogoliubov coefficients defined by 

< = Sa° k ut S^ = a k a° k ut + p k b°£ , (B12) 

We use the relations Eq. (B5) and Eq. (B6) and get 

a k = e w *ch(r k ) , (B13a) 
p k = e - iek e i2(t>k sh(r k ) . (B13b) 

Then one can check that inserting Eq. (B13) in Eq. (18) one recovers Eq. (Bll). 
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